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Many spacecraft, such as the Global Precipitation Mission (GPM) spacecraft, use 


propellant management devices (PMDs) to control liquid propellant within propellant tanks. 
Surface Evolver, an energy minimization routine, is used to find liquid-gas interface shapes 
due to surface tension and propellant center of mass locations for varying initial liquid 
positions, contact angles, and fill fractions for the GPM tank with a PMD. This study shows 
that Surface Evolver can accurately model the liquid-gas interface and liquid center of mass 
location for a complicated tank and PMD. The Surface Evolver results show that the initial 
position of the propellant, the contact angle between the propellant and solid surfaces, and 
the fill fraction of the propellant can drastically affect the liquid-gas interface in the tank 
and thus the center of mass of the propellant. 


cm, 
cm, 

cm, 

dA 

AS Abonar 
As Bboitdy 
ASQionds 
ASQtopcy 
E 


Er 


A> AIS 


g 


Sarsibad 
areqbottfunc 
areciop 


areaiopfunc 


lou voll 


Nomenclature 


lower integration limit along revolve axis 

dry surface area 

isurface area if liquid-gas interface, also referred to as the free surface 
total surface area of tank 

wet surface area 

upper integration limit along revolve axis 


= focal circumference of the circle which is perpendicular to the z-axis 
= x-component of propellant center of mass 


~ y-component of propellant center of mass 


z-component of propellant center of mass 


= differential area of facet 
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differential surface area equation in the x-direction for bottom of sphere 
differential surface area equation in the y-direction for bottom of sphere 
differential surface area equation in the x-direction for top of sphere 
differential surface area equation in the y-direction for top of sphere 
potential energy due to gravity for the propellant 

total energy due to surface tension 

liquid-gas surface tension 

liquid-solid surface tension 


= solid-gas surface tension 


= liquid climb height along the tank wall from center of tank (origin) 
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unit vector in x-direction 

unit vector in y-direction 

unit vector in z-direction 

density of propellant 

surface area of revolved surface 

surface area of bottom half of sphere 

surface area of wetted section of bottom part of the sphere as a function of # 
surface area of top half of sphere 


= sutface area of wetted section of top part of the sphere as a function of A 
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contact angle as measured from the solid to liquid surfaces through the liquid volume 
volume of propellant 

coordinate in the x-direction 

coordinate in the y-direction 

coordinate in the z-direction 
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I. Introduction 


NOWING the center of mass of the liquid propellant is important for attitude control, especially for 
pacecraft which have a large wet/dry mass ratio and/or spacecraft that have large thrusters that can exert large 
moments on the spacecraft. Knowledge of the liquid/gas interface is important for identifying if the outlet areas of 
the tank will be wetted after maneuvers and when using thermal capacitance models to gauge propellant during 
flight. On earth, gravity causes a flat liquid/gas interface to be created in a propellant tank. The liquid/gas interface 
formed in microgravity is caused by surface tension. This results in unintuitive liquid/gas interface shapes, 
especially when a vane-type propellant management device (PMD) is used to orient the propellant at the exit of the 
tank. 

Surface Evolver’ is an energy minimization routine that can be used to model the steady-state liquid-gas 
interface and liquid center of mass position within propellant tanks.*”’ The energies for propellant tanks can include 
surface tension energy and potential energy (body accelerations). Surface Evolver only models the liquid surface, 
and hence converges quickly in most cases. Surface Evolver was developed by Ken Brakke, circa 1991, as part of 
the Geometry Supercomputing Project, which was funded by different academic, government, and private entities." 
It is provided free of charge, with adequate documentation provided by Ken Brakke. Several good review papers 
exist that summarize the work done with Surface Evolver as applicable to propellant tank and PMD design.*"” 

This paper will present methods and results for calculating the liquid-gas interface and liquid center of mass for 
the propellant in the Global Precipitation Mission (GPM) tank with a slightly simplified version of the GPM 
propellant management device (PMD). The liquid initial position in the tank, contact angle between the liquid and 
solid surfaces, and the liquid fill fraction were varied to give a complete view of how the liquid will behave after 
maneuvers. 


Hf. Background 


A. Geometry and Terms 

Figure 1 shows the PMD with its various parts called out. Figure 2 shows the GPM tank and PMD with the axes 
used for this paper, The 
blue lines outline the 






tank wall and the red «JPOP Feservolr tf 

lines outline the PMD. vanes 2 Outlet vanes 
The following terms 

are defined as used in Central vanes 

this paper: contact and 

angle, free surface, communication 

liquid initial position, tube 

Bord number, body 

acceleration, and 

directions associated Bottom reservoir 

with axes. wal 


vanes, sump, and 


The contact angle is 
outlet 


the angle formed 
between a liquid and 
solid. When this angle Figure 1. GPM PMD with different parts labeled. 
is large, the solid 
surface is considered to be non-wetting. When this angle is small, the solid surface is considered to be wetting. A 
zero degree contact angle is considered to be perfectly wetting. 

The free surface is the interface between the liquid and gas in the tank. 

The liquid initial position is defined as the location of the liquid when a simulation is begun. Physically, this is 
the liquid position once a spacecraft maneuver has ended and all sloshing has stopped (or nearly so). 
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The Bond number 
(Bo) is a non-dimensional 
term used to compare the 
relative effect of surface 
tension and a _ body 
acceleration. If the Bond 
number is much smaller 
than 1, the surface tension 
is the dominant fluid 
force. If the Bond number 
is much greater than 1, the 
body acceleration 
becomes the dominant 
force. The body 
acceleration on earth is 
gravity, but in space it is 
the acceleration caused by 
the spacecraft accelerating 
or rotating in a specific Me a / 
direction. kK STE ee 

The directions used for : Mae 
this paper are up, down, 
left, right, forward, and Figure 2. GPM Tank and PMD showing tank dimensions and axes used in 
back. Up is in the positive this paper. 
z-direction, right is in the 
positive y-direction, and forward is in the positive x-direction, with the opposites of these terms being in the opposite 
direction along the same axis. 


39.32 inches 


45 inch 
—— | Tank domes 





B. Low-Gravity Fluid Theory 

The absence of significant body accelerations during much of the spacecraft flight results in surface tension 
forces dominating fluid flow (low Bond number). The lack of body accelerations allows the propellant to be 
modeled by only looking at surface tension. Because surface tension dominates the fluid behavior in this problem, a 
brief review of the equations for calculating the energy associated with surface tension is beneficial to understand 
the results presented in this paper. 

Surface tension is an interaction between different states of matter. There is a liquid-solid surface tension, ;,, a 
solid-gas surface tension, y,,, and a liquid-gas surface tension, Yig- Equation (1)"* shows how the different types of 
surface tension relate in a physical system when the contact angle, @, is known. 


Vsg ~ Vis = Vig Cos @ (1) 


Surface tension is measured in N/m in SI units and can be considered an energy density. Equation (2) shows 
how energy per surface area J/m’ is equivalent to N/m where J is Joules, m is meters, and N is newtons. It is because 
of this relation that sometimes the various surface tensions are referred to as surface energies; however, for this 
paper the term surface tension will be used. 


J _ kg-m? 


m2 s2-m2 





~ (2) 


The energy in a propellant tank problem due to surface tension, E;, can be represented by multiplying the surface 
tensions in Eq. (1) by the surface area with which each surface tension is associated. The different surface areas are 
dry surface area, A,,, wetted surface area, A,,,, and free surface area, Ay. Equation (3) shows the final the result. 


Ey = Aags¥sg + Aws¥its + ArsVig (3) 
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By noting that the total surface area, Az, of the tank is a constant, we can replace either the dry surface area or the 
wetted surface area in Eq. (3) with a term involving the total surface area. The relation between the surface areas is 
shown in Eq. (4). 


Ays = Ar — Aas (4) 
By substituting Eq. (1) and Eq. (4) into Eq. (3) and rearranging terms, Eq. (5) is obtained. 
Er = Aas¥ig COS@ + Ags¥Vig + ArYis (5) 


The surface tensions and total surface area of the tank are constants, so the term, A7y,, in this problem can be 

ignored when doing an energy minimization. Equation (5) is the equation used in Surface Evolver when modeling a 

zero degree contact angle. This is because A,, gives the same contribution as A,, due to the cosine term equaling one. 
By rearranging terms in Eq. (4), and then plugging the result into Eq. (3), Eq. (6)-is obtained, 


Er = —AvwsYia cos @ + AgsVia a ArYsa (6) 


As with Eq. (5), the surface tensions and total area are constant and so the term, A774 is ignored when Equation (6) 
is used in to find the energy minimum in Surface Evolver. 


C. Method for Calculating Surface Energies 
Surface Evolver automatically calculates the free surface area, and with some user programming, can be made to 

calculate the wetted or dry surface area. | 
In the case of a zero degree 

contact angle, the dry surface 

area has the same energy density 

as the free surface area and so 

there is no need to program in the 

dry surface area. Instead, the gas 

bubble is modeled. When the 

bubble touches, or "dries out" the 

tank wall, Surface Evolver 

calculates the surface area of the 

bubble facet, and therefore, the 

dry surface area. All surfaces 


which are not touched by the air 
bubble are considered wetted. Polylines 
Since Eq. (5) does not explicitly : 
contain this quantity it is not Axis about 
important for the calculation of which the 
the total energy. Polylines are 

‘The wetted or dry surface rotated. 


area needs to be calculated 
explicitly within Surface Evolver 
for non-zero degree contact 
angles. The wetted surface area is 
calculated for the cases presented 
in this paper since Eq. (6) is used 
to calculate the total energy. For Figure 3. Polylines that make up the GPM tank outline and axis 
simple shapes, such as spheres, about which the polylines are rotated to create GPM tank surface. 
cylinders, and hemispheres, the 

calculation of the wetted area through a line integral along the three phase boundary is straight forward. The 
example of a sphere can be found in the Surface Evolver manual. 
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GPM uses a cylindrical tank with isotensoid domes. The shape of the isotensoid domes is described within 
Surface Evolver by seven curves that make up a polyline. This polyline is then revolved around the axis that runs 
down the center of the tank. The polylines and rotation axis are shown in Figure 3. 

To illustrate the method used to calculate the wetted energy of the tank wall, and to keep the illustration 
reasonably simple, the example of a sphere is used. Even though the sphere’s wetted surface area can be obtained as 
shown in the Surface Evolver manual, it provides a good example of the method used to find the wetted surface area 
m the GPM tank due to the simplicity of its shape. The sphere is shown in Figure 4 below with the origin of the 
coordinate system at the center of the tank. The sphere has a radius of 1. 






x 
Curve that describes & 
wall outline 
Axis about 
which the 
curve is 
revolved 
a) b) 


Figure 4. Drawing showing the a) curve which is revolved and the b) resulting sphere, 


Step I involves calculating the surface area of a revolve of the polyline that describes the outline of the tank 
wall. Note: this step can be ignored if an equation already exists for the surface area of the tank shape, such as with a 
sphere, but for demonstration purposes we will show the step here. Calculating the surface area, Sareq, Of a revolve 
can be done through Eq. (7)'* where a and b are the lower and upper integration limits along the revolve axis, 


respectively. 
b dx 
Sarea = f, 2mx I + (F)?dz (7) 


Equations (8) and (9) are the result when combining Equation (7) and the equation for the outline of the arc 
shown in Fig. 4. Spreahow and Sareatop are the surface areas for the bottom half of the sphere and the top half of the 
sphere, respectively. When the top and bottom areas are summed the known surface area of a sphere with a radius of 
1, 4x, is found. 








Sareabott = f_,2nVI— 2? |1 + (G5)'dz = 2m (8) 
Sareatop = Jy 2nVI =z? f1 + (dz = 2n 0) 
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Step 2 requires the defining of an equation that gives the local circumference, cf, of the tank as a function of 
distance along the rotation axis. For the purposes of this study the distance along the rotation axis, h, corresponds to 
the height to which the liquid climbs up the wall. Equation (10) shows the local circumference equation for a sphere. 


cf = 2nJ/(1—h)(1 +A) (10) 


Step 3 requires an equation that allows the calculation of the surface area given different liquid climb heights 
along the tank wall. This step involves modifying Eq. (8) and (9) so that the definite integral involves the extra 
variable, 4. Equations (11) and (12) are the resulting equations where Sareaborfunc ANA Soreatopfune are the equations used 
when the liquid climb height is below the z=0 plane or above z-0 plane, respectively. 








h 2 

Sareabottfunc = i 2mvi—z? j1+ ("dz = 2nh+ 2K (11) 
A — 

Sareatopfunc = J 2mvi— 24 j1+ (Faz = 2mh (12) 


Both equations are valid at A = 0. It is important to note that if there are two separate pools of liquid, then the 
direction from which A is 
currently measured (bottom to 
top) would need to be switched 
(top to bottom) for the pool at the 
top of the tank, and hence, each 
pool would have a unique bottom 
and top. 

Step 4 involves combining 
what has been done in the 
previous three steps to create 
equations that Surface Evolver 
can integrate. Figure 5 shows the 
differential surface area, Asa, as a 
function of differential arc length, 
Acf. Equations (13), (14), (15), 
and (16) show how all the 
information above will be put 
together. These equations give the 
differential area, dsa, for the Figure 5. Figure showing Asa and Acf for the sphere. 
bottom and top integrated in the x | 
and y directions. As@pond is for the bottom part of the sphere in the x-direction, ASQiepa. 18 for the top part of the 
sphere in the x-direction, Asazouomdy i8 for the bottom part of the sphere in the y-direction, and AsSQtopa, is for the top 
part of the sphere in the y-direction. 


ct z we 





AsQpoteax = Gam VU — A) +h) Sareabottom) (13) 
Asdtopax = Gap VO — bh) +h) * (Sareatoprunc + Sareabott) (14) 
ASQyottay = aay PGs) ] (Sareabattfunc) (15) 
Asdtopay = Gay VC —h)G +h) = Sareatopfune + Sareabott (16) 


Surface Evolver will solve the equations for each edge that makes up the curve that marks the three phase 
boundary. It will then add up the results to give the final wetted surface area. In practice the surface area is a means 


6 
American Institute of Aeronautics and Astronautics 


of obtaining the surface energy. To convert Eq. (13), (14), (15), and (16) into energy equations, multiply them by the 
ccsine of the contact angle in radians as explained in Section B. 


D. Method for Calculating the Center of Mass 

For spacecraft maneuvers, it is essential to know the center of mass of the propellant within a propellant tank. 
Surface Evolver can easily compute the center of mass of the propellant. This paper. will not present a 
mathematically rigorous derivation of the method used, but will present the author’s methodology for determining 
the correct way to compute the center of mass within Surface Evolver. 

Starting with the equation used to calculate the potential energy due to gravity within Surface Evolver, Eg. (17)’, 
it can be shown how to program Surface Evolver to find the center of mass location of the propellant. Equation (17) 
shows the equation for calculating gravitational potential energy, E, over a facet. The unit vector in the z-direction, 
k, differential area, dA, gravity acceleration, G, and density, p, are included in Eq: (17). Equation (18) shows a 
simplified version of Eq. 17 that includes the z-component of the center of mass, cm,, and the volume, V, of the 
propellant. 


E = G(p) [f=2?k-dA (17) 


E = G(p)(V)(cmz) (18) 


A comparison of Eq. (17) and Eq. (18). shows that the propellant volume is included in the integral in Eq. (17). 
By dividing out V from the integral in Eq. (17), an expression can be obtained for cm,. It should be noted that V is a 
user inputted value in Surface Evolver. Equation (19) is the resulting equation. 


em, = (ff 527k dA)/v (19) 


The process shown above is repeated to obtain the x-component, cm,, and the y-component, cm,, Of the center of 
mass. The resulting equations are Eq. (20) and Eq. (21) with the unit vectors in the x-direction, t, and y-direction, J. 


cm, = (ff $x7i- dA) /V (20) 
em, = (ff 57 dA/v (21) 


For the above equations to work, the entire liquid surface must be facetted, including the liquid-solid interface. 
The exception to this rule is when a tank wall is a flat surface at the x = 0, y = 0, and/or z = 0 planes for cm,, CMy, 
and/or cm,, respectively. In these cases the walls at the zero plane can remain unfacetted and still give an accurate 
value for the center of mass, however, the rest of the walls should be facetted. The simulation will run faster by 
keeping the number of faces to a minimum, so keeping as many walls unfacetted as possible is recommended. 


If. Results an Discussion 


A. Zero Degree Contact Angle Results 

The first set of results assumes design conditions for the GPM tank at launch. GPM’s tank was designed to have 
a zero degree contact angle with hydrazine. The approximate fill fraction at launch is 70%. The initial position of the 
liquid propellant is varied to show how the hydrazine propellant would behave if the launch pushed the liquid to 
different regions of the tank. The initial propellant positions assumed that the PMD would have filled up, as per its 
design.., 

Figure 6 shows the final solution for the zero degree contact angle cases. The white bubbles represent the gas 
bubbles. The liquid propellant is wherever the bubble is not. 
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The imitial position results shown in Fig. 6 give some interesting insight into the fluid behavior. First, the bubble 
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Figure 6. Liquid-gas interface results for zero degree contact angle 
cases at a fill fraction of approximately 70% with varying initial 


positions (see label above each subfisure for the initial position). 


where the bubble touches 
the wall, gives a solution 
where part of the wall is 
dried out. Fourth, the 
total energy is similar for 
cases where the bubble is 


will not automatically go back to the 
center of the tank. Instead, the 
bubble tends to remain near its 
initial position. Second, the final 
resulting bubble position does not 
dry out the PMD vanes, even if its 
initial position has it drying out 
some of the vane. Third, if the 
bubble touches the wall of the tank 
it will dry out a part of the tank wall. 
This third observation goes against 
the assumption that a perfectly 
wetting liquid will wet all the solid 
surfaces. The explanation for this 
can be found in Eq. (5). Equation (5) 
shows that it is not only the dry 
surface area, but also the free 
surface area, that contributes to the 
total energy in the problem. The 
combination of dry surface area and 
free surface area, found in cases 


Table 1. Center of mass (CM) resuits in inches (in) for zero degree contact 
angle cases at a fill fraction of approximately 70% with varying initial positions. 


initial Bubble Pasition | Xem{in) | 






completely Eee completely Centered Q 0 
around the center post, Top-Center 0 0 
and where the gas bubble cree -0.03 0 
does not completely wrap ites 0 -2 48 
oer ease = eect 2.48 0.05 
seco € 0 
ee ‘s more Tee Peeks a ele z02 
because any significant Center-Fromt 0 
acceleration would push [MAGUS TGieiians 0.16 -3.75 
the bubble away from the Half-Right-Off-center -1.17 -3.58 
center post. Half-Top-Right -1.56 -3.36 
sis Eee end ee am 383 Hue 
y Peace area 3.64 0.98 


we see in Fig. 6. As 
expected, the center of 


an a) 





Ee cM Magnitude (in) 


0.01 0,01 
-2.07 2.07 
2.17 2.17 
-1.81 3.07 
1.73 3.02 

0 2.63 
0.01 2.69 
-0.09 3.75 
0.02 3.77 
-0.09 3.82 
-0.03 3.77 
-0.40 3.79 


mass is offset more from the center of the tank if the bubble is not wrapped around the center post. 


B. Nou-zero Degree Contact Angle with PMD 

The third set of results examines the liquid propellant behavior in the GPM tank and PMD for varying fill 
fractions, contact angles, and initial positions. Throughout spacecraft operation, the fill fraction in the tank decreases 
and the position of the propellant changes. The contact angle changes when contamination enters the tank or 
prepellant. Figure 7 shows the results with the liquid initial position being almost equal on the top and bottom of the 
tank. The fill fractions are varied between 3% and 70% and the contact angle is varied between 10 deg and 60 deg. 
The dark area represents the liquid position. The faces that break up the liquid-gas interface are transparent and so 
the darker areas within the shaded area show where the air bubble goes below the top most part of the liquid within 
the top or bottom liquid pools. 
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The results 
show two 
interesting 


trends. First, as 
the contact angle 
increases, the 
wetted surface 
area on both the 
vanes and tank 
wall decreases. 
Second, as the 
fill fraction 
increases, the 
region in which 
three-phase 
contact line 
occurs changes. 
These behaviors 
change the shape 
of the liquid gas 
interface and 
hence the 
location of the 
center of mass. 


more _ liquid 
propellant at 
the top of the 
tank, the 
ability of the 
PMD to move 
the liquid to 


the outlet of the tank can be studied. Figure 8 shows the results 
initialized at the top of the tank, and the contact angl 
Figure 8 shows that as the contact an 
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oo 
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Figure 7. GPM tank and PMD liquid-gas 
interface case results with fluid equally 
distributed between top and bottom of tank, 
varying contact angles, where a) 10 deg, b) 30 
deg, c) 60 deg, d) 10 deg, e) 30 deg, 60 deg, 2) 
10 deg, h) 30 deg, i) 60 deg, j) 10 deg, k) 30 deg, 
I) 60 deg, m) 10 deg, n) 30 deg, and o) 60 deg, 
and varying fill fractions (see legend). 
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Figure 8. GPM Tank and PMD 
liquid-gas interface case results at a 
fill fraction of 20% and varying 
contact angles of a) 10 deg, b) 30 
deg, c) 60 deg with more propellant 
initialized toward the top of the 
tank than at the bottom of the tank. 


Table 2. Center of mass (CM) results for non- zero degree contact angle cases at a fill 
fraction of approximately 20% with varying initial positions. 






Pn CM Magnitude (in) 


8.00 1.35 
1Q.59 6.90 
-18.06 L161 


at a fill fraction of 20%, with more propellant 
e varied from 10 deg to 60 deg. 
gle increases, the air bubble results are less centered about the center of the 


tank. Table 2 shows numerically these differences in center of mass location. - 


Understanding how liquid propellant behaves in a low- 
propellant tanks. Surface tension forces cause the liquid prop 


IV. Conclusions 


cenier of mass and the liquid-gas interface shape. 


Surface Evolver can be used to model the li 
liquid propellant. The results in this report 


gravity environment is important for the design of 
ellant to orient in unintuitive ways, thus changing the 


quid-gas interface shape and hence the center of mass location of the 
show that a complicated PMD and non-standard tank shaped tank, such 


as GPM’s PMD and propellant tank, can be modeled for both non-zero and zero contact angle cases using Surface 
Evolver. The liquid-gas interface and hence the center of mass of the liquid propellant is affected by the liquid fill 
fraction, position of the liquid when thrusters (or other acceleration) is turned off, and the contact angle between the 
liquid propellant and the solid walls and vanes. As the contact angle decreases the liquid will wet more of the solid 
surfaces in the tank, but may not wet all the solid surface even at zero degree contact angle. 


Future work in this area will include modifying 
face initializing techniques. This will enable higher 


Surface Evolver’s wetted energy models and developing better 
fidelity results and more stable simulations. 


Physical experiments are essential to showing that the Surface Evolver models are giving correct results. The 


microgravity experiments that could validate Surface Evolver inc 


drop-towers, with sub-orbital flights, or on-orbit flights. 
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